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Abstract The origin of the regular satellites ties directly to planetary formation in that the satellites form in 
gas and dust disks around the giant planets and may be viewed as mini-solar systems, involving a number 
of closely related underlying physical processes. The regular satellites of Jupiter and Saturn share a number 
of remarkable similarities that taken together make a compelling case for a deep-seated order and structure 
governing their origin. Furthermore, the similarities in the mass ratio of the largest satellites to their pri- 
maries, the specific angular momenta, and the bulk compositions of the two satellite systems are significant 
and in need of explanation. Yet, the differences are also striking. We advance a common framework for the 
origin of the regular satellites of Jupiter and Saturn and discuss the accretion of satellites in gaseous, circum- 
planetary disks. Following giant planet formation, planetesimals in the planet's feeding zone undergo a brief 
period of intense coUisional grinding. Mass delivery to the circumplanetary disk via ablation of planetesimal 
fragments has implications for a host of satellite observations, tying the history of planetesimals to that of 
satellitesimals and ultimately that of the satellites themselves. By contrast, irregular satellites are objects 
captured during the final stages of planetary formation or the early evolution of the Solar System; their dis- 
tinct origin is reflected in their physical properties, which has implications for the subsequent evolution of 
the satellites systems. 
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The satellites of Jupiter and Saturn have received renewed interest as a consequence of the Cassini mission to 
Saturn, which has resulted in a growing body of literature concerning both irregular (see Jewitt & Haghigh- 
ipour (2007) for a review) and regular satellites (see Estrada et al. (2009) for a review; and Canup & Ward 
(2009) for a different take). We make no claim to completeness - the aim of this submission is not to provide 
a comprehensive review of the field, which has already been attempted in those publications, but to describe 
key physical processes and observational constraints. While there is unambiguous observational evidence 
that the regular and irregular satellites are distinct populations, there exist at least indirect ties between them 
when it comes to the processes involved in the capture of objects from heliocentric to circumplanetary orbits. 
There is also a hint of a direct connection in terms of the locations of the innermost irregular satellites and 
the size of the disk formed by gas inflow through a gap (^ Rh/5', Mosqueira & Estrada (2003a); Estrada et 
al. (2009); Ayliffe & Bate (2009)); irregular satellites captured into orbits with smaller semi-major axes may 
be removed by gas drag from the outer disk and possibly contribute mass to the regular satellites, though this 
link remains to be quantitatively explored. 

It is also worth noting that the processes shaping the formation of satellite systems are strongly coupled to 
those governing planetary formation and the evolution of the early Solar System. In fact, planetary forma- 
tion sets the initial conditions for satellite formation which, therefore, cannot be treated in isolation. The 
forming satellite systems can exchange material with the outer Solar System, receiving solids mainly in the 
form of planetesimals captured or ablated due to the interaction with the circumplanetary gas disks, and 
even scattering satellitesimals out into heliocentric orbits. Moreover, the regular satellites, once formed, are 
still subject to exchanges with the outer Solar System, as dramatically illustrated by the disruptive capture 
of Triton in Neptune's satellite system. Below we discuss the formation and properties of both regular and 
irregular satelUtes in the framework of such exchange processes. 



1 Formation of the regular satellites of giant planets 

In the core accretion model of planet formation (Bodenheimer & Pollack, 1986) a core must first form by 
accretion of planetesimals. In this mode of planet formation most of the mass of solids is taken to reside in 
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planetesimals of size ~ 10 km (Wetherill & Stewart, 1993). Planetesimals may dissolve in the envelopes of 
the forming giant planets, thus enhancing the planet's metallicity (Pollack et al., 1986). Most of this high-Z 
mass deUvery takes place before the cross-over time, when the mass of the gaseous envelope grows larger 
than that of the core. This stage is then followed by a dilution during runaway accretion, depending on the 
amount of gas accreted during this brief phase of the planet's growth (Pollack et al., 1996). As the giant 
planet reaches its final mass, planetesimals in its feeding zone undergo coUisional grinding. In the Jupiter- 
Saturn region, the colUsional timescale for kilometer- sized objects is similar to the ejection timescale ~ lO"* 
yrs, so that a fraction of the mass of solids will be fragmented into objects smaller than 1 km (Stern & 
Weissman, 2001; Charnoz & Morbidelli, 2003). The coUisional cascade facilitates planetesimal delivery to 
the circumplanetary disk (in ~ 10^ yrs) because smaller planetesimals are easier to capture. 



1.1 Satellite migration and growth 

Mosqueira & Estrada (2003a) start with a disk populated by satellite embryos (^ 10^ km) and satellitesimals. 
These authors adopt a bimodal size distribution or two-groups method, as has been done in a number of 
pubhcations that treat protoplanetary growth by the binary accretion of planetesimals (e.g., Pollack et al. 
(1996); Goldreich, Lithwick & Sari (2004)). There are close parallels, but some differences as well. In both 
the planet and satellite cases, the migration of the embryos is dominated by the tidal interaction with the 
gas disk, whereas planetesimals and satellitesimals migrate due to aerodynamic gas drag. But in the satelUte 
case the transition from one type of migration to the other overlaps and leaves no size unaffected, i.e., 
all object sizes from kilometers to full-sized satellites undergo migration. Another difference involves the 
capture cross-section of migrating planetesimals versus satellitesimals. Satelhte embryos have a larger ratio 
of physical radius to Hill radius (^ 0.1 at the location of Ganymede) than do protoplanets (~ 10~^ at the 
location of Jupiter). As a result, the impact probability of passing objects is larger (cf. fig. 9 from Kary, 
Lissauer & Greenzweig (1993)). Also, the resulting gravitational focusing factor {Fg ~ 10) for satellites is 
never as large as it is for planet embryos during their runaway phase, before they grow large enough to "heat" 
up their own food and their growth switches to oligarchic accretion^ 

1.1.1 Timescale for disk formation and cooling 

While the gas accretion rate is high the planet's envelope takes up a fraction of the planet's Roche lobe. 

Thus, there is agreement in the literature that the formation of icy satellites close-in to the planet must await 
for gas accretion onto the planet to wane (Coradini et al. (1989); Makalkin et al. (1999); Ayliffe & Bate 
(2009)). The timescale for Jupiter and Saturn to clear a gap in between (t ~ 10"* years or less in a compact 
initial configuration) sets the timescale for the end of gas accretion^. The timescale for envelope collapse is 
given by the Kelvin-Helmholtz time of ~ 10^ years (Hubickyj, Bodenheimer & Lissauer, 2005), following 
which a circumplanetary disk will form. We stress that even after envelope collapse continued inflow onto 
the disk may keep it too hot for the concurrent accretion of close-in, icy satellites (see Estrada et al. (2009) 
and Ayliffe & Bate (2009) for further details). Once the circumplanetary disk forms and gas inflow dies down 

' At this stage, satellitesimals should be interpreted as those objects which determine embryo growth, containing a significant 
fraction of the mass of solids in the disk (though other sizes may also contribute to the mass distribution). 
^ Provided the protoplanetary disk at the location of Jupiter and Saturn is weakly turbulent (Bryden et al., 1999). In that case 
the gas inflow through the gap ceases. 
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(as the giant planets open a combined gap), the timescale for the disk to cool sufficiently for ice to condense 
depends on the disk opacity. The Rosseland mean thermal opacity for the gas is Kgas ^ 10^^ cm^g^\ 
while that of micron-sized dust grains can be of order ~ 1 cm^ g^' (Lunine & Stevenson, 1982; Chiang 
et al., 2001). As particles grow, the opacity of the nebula decreases. A sustained high opacity may keep 
the subnebula weakly turbulent and hotter Conversely, rapid coagulation allows the subnebula to become 
optically thin, cooling quickly and allowing turbulence to subside. Indeed, while the incremental growth of 
grains and dust may take place irrespective of the level of nebula turbulence, growth past the decoupling size 
likely requires very low levels of turbulence (a < 10^^) (Youdin &. Shu, 2002; Cuzzi & Weidenschilling, 
2006), which allows for a cool subnebula where condensation of volatiles takes place (Mosqueira & Estrada, 
2003a). Following envelope collapse, ablation of meter-to-kilometer planetesimal disk crossers (see Estrada 
et al. (2009) and references therein) and other mechanisms enhanced the concentration of solids in the 
satellite disk from a starting condition with condensable content enhanced by a factor of 3 — 4 from cosmic 
mixtures, as suggested by the high-Z content of giant planet envelopes (Atreya et al., 1999). Note that the 
circumplanetary disk and the giant planet envelope are both enriched in solids as a result of a common 
mechanism, i.e., planetesimal ablation (see Fig. 1 upper left panel; Mosqueira, Estrada & Charnoz (2009)); 
however, the quantitative degree of enrichment may differ substantially between the two^. Given a density 
of ^ 1 g cm^^, a planetesimal of size < 1 km encounters a gas column equal to its mass when crossing the 
circumplanetary gas disk. Such a population of objects is hkely to result from the fragmentation cascade 
following giant planet formation (e.g., Charnoz & MorbidelU (2003)). Thus, we consider a soUds-enhanced, 
quiescent disk with optical depth T > 1 and gas surface density ^10"* g/cm^. The disk cooling time is then 
given by Xkh ~ ^gCTc/ {IosbT^) ~ lO^ilSOK/Tcyx yrs, where Osb is the Stephan-Boltzmann constant, C 
is the specific heat, and Tc and are the midplane and photospheric temperatures. 

1.1.2 Disk scales 

Gas flowing into the protoplanet's Hill sphere forms a circumplanetary disk. One can obtain an estimate 
of the expected disk size by assuming that in-falling gas elements conserve specific angular momentum. 
Assuming that prior to gap-opening the giant planet accretes gas with semi-major axis originating within 
^Rh from its location, centrifugal balance yields a characteristic disk size of ~^///50 (Cassen & Pettibone, 
1976; Safronov et al., 1986; Mosqueira & Estrada, 2003a). For Jupiter and Saturn these radii are located 
close to the positions of Ganymede and Titan, respectively. On the other hand, after gap-opening accretion 
continues through the planetary Lagrange points. In this case, the estimated characteristic disk size formed 
by the inflow is significantly larger, of the order of ~ Rn/'i (Mosqueira & Estrada (2003a), but note that 
this paper settles on an outer disk size of ^ Rh/^, based on the locations of irregular satellites). In the 
absence of a gap, 3D effects allow low angular momentum gas to be accreted directly onto the planet, or in 
a compact disk of characteristic size < O.IRh (D'Angelo, Kley & Henning, 2003), whereas in the presence 
of a well-formed gap a 2D treatment applies and the resulting disk size will be > O.IRh- D'Angelo, Kley 
& Henning (2003) conclude that, for nominal disk parameters, two dimensional computations are reliable 
for secondary to primary mass ratio ju > lO""*, such that the Hill radius of the protoplanet is larger than the 
disk scale-height. Estrada et al. (2009) and Ayliffe & Bate (2009) show that the specific angular momentum 
of the inflow through the gap is about a factor of 3 — 4 larger than that of the Galilean satellites 1.1 x 10^^ 
cm^/s, which means that the gas will achieve centrifugal balance at a radial location of ^ 100/?/. This disk 



^ The runaway gas accretion phase is expected to result in a dilution of the solid content of the planet's envelope, as most of the 
mass of solids is in planetesimals that are not coupled to the gas (see Estrada et al. (2009) for a detailed discussion). 
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is compact compared to the Hill radius, but it is extended in terms of the locations of the regular satellites, 
and possibly linked to the location of the irregular satellites. 

1.1.3 Model parameters 

Mosqueira & Estrada (2003a) divide the circumplanetary disk into inner and outer regions. For Jupiter, they 
compute the solids-enhanced minimum mass (SEMM) gas densities in the inner and outer disks based on the 
sohd mass required to form lo, Europa (both re-constituted for missing volatiles) and Ganymede in the inner 
disk and Calhsto in the outer disk. Inside of the centrifugal radius the surface gas density Zg is 10^ — 10^ 
g cm^^ (corresponding to pressures ^0.1 bar), which yields an optical depth due to absorption by hydro- 
gen molecules (assuming dust coagulation) Ty ~ ZKgas ~ 1. Outside the centrifugal radius ~ 30/?/, the 
gas surface densities are in the range 10^ — 10^ g cm~-^, which results in a low optical depth in the range 
T,, ^ 0.01 — 0.1. These authors apply the same procedure to the circumplanetary disk of Saturn by employ- 
ing the masses of Titan and lapetus to set the inner and outer disk masses, respectively, and choose a simple 
model where the gas density follows a simple I /a dependence inside of ~ 20Rj and outside of ~ 26Rj. 
The transition region has a width of (at least) 2Hc, where is the subnebula scale-height at the centrifu- 
gal radius. This choice ensures that the gradient in gas density is not so steep as to lead to a Rayleigh-Taylor 
instability (e.g. Lin & Papaloizou (1993)). The use of a two component subnebula hinges on the assumption 
of decaying turbulence as Roche-lobe gas inflow ebbs. But it should be noted that detailed simulations of 
the formation of such a disk remain to be carried out'*. We tie the disk's temperature profile to the planetary 
luminosity at the tail end of giant planet formation (Hubickyj, Bodenheimer & Lissauer, 2005). For now, we 
use a heuristic temperature profile of the form T = 3600Rj/a for Jupiter (e.g., Lunine & Stevenson (1982)) 
and T = 2QQQRs/a for Saturn. The outer disks of Jupiter and Saturn have roughly constant temperatures in 
the range of 70 — 130 K for Jupiter and 40 — 90 K for Saturn, depending on solar nebula parameters. 
Inner disk: In analogy to planetary accretion, we begin by calculating characteristic object sizes in such a 
disk. The characteristic size over which the self-gravity of the solids overcomes the shear is /i = r/j, where 
/i = (mi/nZsy^^, Eg ~ lO'^ g/cm^ is the surface density of sohds, and = a(mi/3M/>)*/^ is the Hill radius 
of a satellitesimal of mass nii. The second characteristic mass can be obtained from the condition I2 = rn, 
where h = mi/i^'JtaLs) is the distance over which the mass rm in a disk annulus with semimajor axis a 
can force close encounters with passing satellitesimals, and rn = a(m2/3M/))^/^ is the Hill radius of an 
embryo of mass m2. For object densities of ^ 1 g/cm^, these characteristic masses correspond to sizes of 
ri ~ 1 — 10 km and r2 ~ 1000 km. We now consider a gaseous disk with surface density Eg ~ 10"* g/cm^ and 
aspect ratio ///a ~ 0.1 (or a temperature of ~ 250 K at Ganymede [l5Rj], and ~ 100 K at Titan [20/?,]), and 
obtain a Toomre parameter Q ^ 100, which shows that self-gravity is unimportant. The particle decoupling 
size is £2ts = l6psrpVk/ {3CDEgcQ) ~ 1 ^ ~ 10 m, where ts is the stopping time, is the Keplerian 
velocity, Cd ~ 1 is the drag coefficient^ (see e.g. Adachi, Hayashi & Nakazawa (1976)), and c is the sound 
speed. The gas drag time for a satellitesimal of size ri located at l5Rj is T^rag = 1 /2{vk/c)^ts ^ 100 years. 
The gas drag time for an embryo of size r2 is ~ 10^ years. Such an embryo will also migrate due to the tidal 
interaction with the gas disk in a similar timescale T/ w QQa^c^ / {ZgG^m2) ~ 10^ years, where G is the 
gravitational constant and C/ < 1 in a 3D disk (Tanaka, Takeuchi & Ward, 2002; D' Angelo, Kley & Henning, 
2003; Bate et al., 2003). Thus, the crossover size between the regimes of gas-drag and tidal-torque migration 

'* Here we rely mostly on the overall gas surface density of the subnebula, which we justify on observational grounds. 
^ There is an issue whether this drag coefficient should include a correction in cases where the satellite can gravitationally 
perturb the gas flow (e.g., Canup & Ward (2002)). Ostriker (1999) argues that this enhanced drag is already incorporated into 
the migration rate due to differential tidal torques. 
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is ~ 500 — 1000 km. The Type I migration timescale for Ganymede is ~ 10"^ years, which is comparable to 

its formation time given by the gas drag time of satellitesimals of size ^100 km. 

The question arises whether the satelhte growth times are faster than the migration times. Ignoring gravi- 
tational focusing, the Safronov binary accretion timescale for embryos is Tacc ~ Psri/i^s^) ~ 10^ years, 
which is faster than either the Type I or gas drag migration for such objects in the SEMM disk. But note 
that if all the satellitesimals resided in objects of size ri ~ 1 km, then the time to clean the disk by gas 
drag migration would be ~ 100 years. The mass growth timescale due to dust and rubble sweep up as a 
result of differential drift is = 4-psr2/ {3ppAvp) ^ Psr2a/{EgQ){a/H) ^ 10^ years, where Avp ^ 77 and 
^^ = {vic — Vg)/vk is the fractional difference between the local Keplerian velocity and the gas velocity 
Vg due to gas pressure support. This estimate assumes that the density of the dust and rubble layer settles 
to a state such that Pp « Pg, where Pg is the nebular gas density, which is consistent with a quiescent disk 
stirred only by weak, local production of turbulence at the shear layer. Thus, satellitesimals grow as they 
drift: We interpret Hyperion (whose radius is 135 ± 4 km, Thomas et al. (2007) ) as a satelUtesimal that was 
captured into resonance by gas drag inward migration (see section 1.2.1). For objects of this (^ 100 km) 
size, satellitesimal migration, and satellite migration and growth all take place on a similar timescale. This 
is because satellite growth is regulated by the rate at which the embryo's feeding zone is replenished by 
satellitesimals or embryos. Also, the Type 1 migration timescale of satellites is comparable to the gas drag 
migration timescales of such satellitesimals (cf. Titan-Hyperion). 

Outer disk: Further out the characteristic sizes are smaller despite the increase in semi-major axis. The rea- 
son is that the surface density is taken to drop outside the centrifugal radius. Typical sizes in the outer disk 
for £^ ~ 10 g/cm^ are ri ~ 1 km and ~ 100 km. The Safronov binary accretion time to form an embryo of 
size ^2 in the outer disk is ~ 10^ — 10^ years. The time for gas drag to clear the outer disk of such embryos 
is also 10^ — 10^ years, which Mosqueira & Estrada (2003a) tie to Callisto's formation time. 

1.1.4 Satellite survival 

So far we have implicitly assumed that the relatively massive gas component of the circumplanetary disk 
is static and unaffected by the evolution of the solids. However, the SEMM model of Mosqueira & Estrada 
(2003b) relies on gap-opening for satellite survival, which means that at some point during their migration 
and growth satellites become sufficiently massive to open a gap in the gas. There are uncertainties involved 
in estimating the value for the gap-opening masses, such as the role of 3D effects and stratification, but they 
are unlikely to greatly affect the inviscid gap-opening criterion of Rafikov (2002) (see also Li et al. (2009)), 
which is based on damping of 2-D acoustic waves by wave steepening (Goodman & Rafikov, 2001). For 
Jupiter and Saturn Rafikov's criterion yields gap-opening masses at l5Rj and 20Rs in disks with /; ~ 0.1 
equal to the masses of Ganymede and Titan for disk surface densities of ~ 2 x lO'* g/cm^ and 1 x 10^ 
g/cm^, respectively. It should be noted that the gap-opening mass depends on the semi-major axis a. Also, 
we assume a surface density profile Zg 1/a, and an extended giant planet following envelope collapse 
with radius Rp ^2Rj (e.g., Hubickyj, Bodenheimer & Lissauer (2005)), where .Ry is Jupiter's radius. As a 
satellite migrates inwards, it finds it easier to open a gap for a given surface density, but the surface density 
is taken to increase closer to the planet. While the timescale for accretion decreases closer-in, the satellite's 
feeding zone becomes smaller, and its inward migration speeds-up. Also, both the migration rate and the 
gap-opening condition are sensitive to the aspect ratio H/a. 
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1.2.1 Capture of Hyperion into a 4:3 resonance with Titan 

The origin of Hyperion in the 4 : 3 mean-motion resonance with Titan presents a significant challenge. A 
tidal origin of resonance capture as may apply to Galilean satellites (Showman & Malhotra, 1997) does not 
apply to the case of Titan and Hyperion. Given Titan's size and distance from Saturn, significant expansion 
of its orbit would require Saturn's dissipation parameter Q to be much lower than the lower hmit set by 
the proximity of Mimas (Lee & Peale, 2000). Work by Weidenschilling & Davis (1985) showed that the 
combination of gas drag and perturbations due to mean motion resonances (MMR's) with a planet can 
have important consequences for the orbital evolution of small planetesimals. Kary, Lissauer & Greenzweig 
(1993) obtained a critical value for the drag parameter at which resonance trapping by a protoplanet breaks 
down. Malhotra (1993) points out that resonance trapping is vulnerable to planetesimal interactions. This 
issue also arises in the satellite case, as the capture of Hyperion in a 4 : 3 MMR with Titan illustrates. 

1.2.2 Callisto's internal state 

A challenging result is that the Galileo mission moment of inertia data are consistent with a fully differ- 
entiated Ganymede, but only a partially differentiated Callisto (Anderson et al., 2001). While it is possible 
that non-hydrostatic effects in Callisto's core could be large enough to allow for complete differentiation of 
Callisto (see e.g. McKinnon (1997); Stevenson et al. (2003)), its internal state may be a result of its accretion 
history. In this view, Callisto forms by accreting the volatile-rich condensables present in the extended, low 
density outer disk; its long formation timescale is tied to the disk clearing time, which is the time it takes for 
gas drag to clear the circumplanetary disk of soUds (Mosqueira et al., 2001; Mosqueira & Estrada, 2003a; 
Ahbert, Mousis & Benz, 2005). That is, accretion in the sparse, extended outer disk takes a sufficiently long 
time that the heat of accretion can escape, and the satelUte may not differentiate fully. Calhsto could then 
be said to be the result of slowly assembling thousands of volatile rich, "cold" embryos. Thus, the issue of 
Callisto's state of differentiation is Unked to the observation of the empty space outside Callisto and inside 
the irregulars. A related question can be posed for Saturn: namely, why is lapetus stranded far from the 
planet? In each case, it is natural to expect that gas drag clearing of swaths of the circumplanetary disk is in- 
volved in the answer^. lapetus' separation from Titan provides strong indirect support for a two-component 
subnebula: a dense iimer region roughly out to the centrifugal radius near Titan, and a lower surface density 
tail extending out to as far as irregular satellite Phoebe. Similarly, Callisto derives its mass from a more 
extended region than does Ganymede, which can explain its partially differentiated state without resorting 
to fine-tuning poorly known parameters; that is, in a two component subnebula it is natural to expect both 
that a) some satelUtes form in the outer disk and others in the inner disk, and b) that satellites forming in the 
outer disk take significantly longer to accrete than those forming in the inner disk. 

The issue arises whether typical impactors would bury heat in Callisto so deeply that it would lead to full 
differentiation of this satellite (McKimion, 2006). Mosqueira & Estrada (2003a) estimate that outer disk 
satellitesimals could grow as large as several hundred kilometers as they migrate inwards due to gas drag. 
These sizes are computed by assuming that all colUsions among satelUtesimals are accretionary. However, 
as satellitesimals approach a satellite's feeding zone they are excited to velocity dispersions well in excess 
of their escape speeds, which is likely to result in non-accretionary, disruptive collisions (as Hyperion's ir- 



* Note that the capture of Hyperion into resonance also involves gas drag inward migration. 
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regular shape may attest to). Thus, it is reasonable to expect that Hyperion-like or smaller sateUitesimals are 
typical impactors during the late-stage accretion of both CaUisto and Titan. 

1.2.3 Compositional constraints 

Regular satelhtes may provide a probe of the compositional and thermal state in the subnebula at the tail end 
of giant planet formation. Early models propose the formation of satellites in the circumplanetary nebulae 
of Jupiter and Saturn (Pollack & Reynolds, 1974; Pollack et al., 1976). These models envision a condensa- 
tion sequence analogous to that of the nebula. Still, the physical conditions in the circumplanetary disks are 
thought to be sufficiently different from those of the nebula that a sharp contrast has been drawn between ob- 
jects forming in the outer nebula and giant planet subnebulae (Prinn & Fegley, 1981). Yet, such ideas retain 
a phenomenological character, as the nature of the interaction between the two environments remains poorly 
understood. The compositional gradient of the Galilean satellites may provide a link to the environment in 
which they form (see Estrada et al. (2009) and references therein). The similarities in the bulk properties of 
the regular satelhtes of Jupiter and Saturn strongly favour a unified framework for their origin; yet, the inner, 
icy satellites of Saturn exhibit no clear compositional trend (and may have been collisionally disrupted^). 
In fact, recent Cassini results make a strong case that the medium-sized Saturnian satellites have densities 
(Jacobson et al., 2006) that preclude a solar composition ice/rock ratio (e.g., Johnson & Lunine (2005) and 
references therein; Mosqueira & Estrada (2005)). On the other hand, a persistently hot subnebula can prevent 
the condensation of volatiles close to Jupiter (Pollack & Reynolds, 1974; Lunine & Stevenson, 1982). Thus, 
the thermal conditions under which coagulation and growth take place may hold a key to the formation of 
the regular satellites. 

Compositional Constraints of the Outer Regular Satellites: Mosqueira, Estrada & Charnoz (2009) focus on 
the large, outer regular satelhtes of each satelhte system: Ganymede and CaUisto in the case of Jupiter, and 
Titan and lapetus for Saturn. For objects the size of lapetus or larger, the porosity is likely to be small not 
only because the internal pressure is large enough to close pore spaces, but also because the presence of 
short-lived radioactive nuchdes heats the interior causing ice to flow. For such large satelhtes, densities can 
be interpreted in terms of rock/ice fractions. lapetus' low density, and correspondingly low rock/ice fraction, 
presents a puzzle when compared to the other three satellites, each of which is roughly 50% ice and rock. In 
turn, the rock/ice fractions for Ganymede, Titan and CaUisto are comparable to that of (captured) Saturnian 
irregular satellite Phoebe (or even smaller if Phoebe were porous). 

Tying the properties of solar nebula planetesimals to subnebula sateUitesimals involves three distinct aspects: 
first, characterising the properties of the first generation of planetesimals in terms of sizes and degree of het- 
erogeneity; second, quantifying the coUisional evolution of the planetesimal swarm following giant planet 
formation; and third, delivering planetesimal fragments to the circumplanetary disk. Planetesimal break-up 
in tandem with delivery via ablation of planetesimal fragments crossing the subdisk provides a framework 
for understanding the mass budget and compositions of regular satellites (compared to that of solar nebula 
planetesimals). In particular, ablation can result in fractionation, and account for the observed density of 
lapetus provided this satelhte formed in situ (Mosqueira & Estrada, 2005). For this to work planetesimals of 
size ^10 km need to be partially differentiated, which indicates that the first generation of planetesimals in 
the Jupiter-Saturn region (and possibly beyond) incorporated significant quantities of ^^Al. 
Thermal ablation occurs because friction of the fast moving heliocentric interloper as it crosses the circum- 

^ Indeed, a number of factors, such as the stochastic compositional gradient of the inner Saturnian satelhtes, the majestic icy 
rings, and Titan's isolation and eccentricity, argue in favour of colhsional processes (see Mosqueira & Estrada (2005) for an 
intriguing, though intricate, possibility). 
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planetary gas disk heats up the body. At low gas densities impactors lose mass and energy through ablation. 
For large kilometer-sized bolides mechanical destruction rather than thermal ablation may dictate the fate of 
the object (see Mosqueira, Estrada & Charnoz (2009) and references therein). Ablation of planetesimal frag- 
ments < 1 km may deliver the bulk of the soUds needed to form the satellites. For objects in that size range, 
the heat transfer coefficient C// ^ 0.1 may be reasonably obtained from observations of terrestrial meteorites 
(Bronshten, 1983). The rate at which energy is transferred to the planetesimal is given by £ ~ c/fp v^, where 
p is the gas density, v is the speed of the planetesimal through the gas, and some planetesimal flattening 
(which increases its cross section) takes place. One can then obtain an estimate of the surface temperature 
Ts by balancing this heating and the radiative cooling. Ablation thus results in delivery of material to the 
circumplanetary disk. Gas drag may also result in capture of material. To quantify these processes we first 
need to characterise the properties of the circumplanetary disk. 

In the SEMM model one expects a factor of ~ 10 enhancement in solids over cosmic mixtures, resulting 
in a gas surface density of ~ 10"^ g cm~^, which is consistent with and quantitatively constrained by the 
gap-opening condition for Ganymede and Titan in a inviscid disk with aspect ratio '--^0.1 (Rafikov, 2002; 
Mosqueira & Estrada, 2003b), and the Type I migration of full-sized satellites in such a disk (Mosqueira & 
Estrada, 2003b; Bate et al., 2003). Following envelope collapse, planetesimal ablation of meter-to-kilometer 
disk crossers (Mosqueira & Estrada, 2005) and other mechanisms enhanced the concentration of solids in 
the satelhte disk. If we use a gas surface density of 10"* — 10^ g cm"^ for the Jovian and Saturnian subnebulae 
then a planetesimal of density ~ 1 g cm~^ will encounter a gas colunrn equal to its mass if its radius is in 
the range 0.1 — 1 km. Following giant planet formation collisional fragmentation may replenish the 0.01 — 1 
km size range, even if most of the mass was originally in objects > 1 km. This means that a significant 
fraction of the mass in planetesimals crossing the gas disk may be ablated. Such a planetesimal may deposit 
a significant fraction of its mass in the gas disk. 

In Figure 1 we show (see Mosqueira, Estrada & Charnoz (2009) for further details) the mass delivery by 
ablation and gas drag capture of planetesimal fragments crossing the SEMM circumplanetary disk. These 
results are obtained assuming a differentiated population of icy and rocky interloper fragments with a rock 
fraction of 30% by mass, and placing ~ 10 Earth masses in between Jupiter and Saturn. The bottom left 
panel shows that it is possible to deliver enough mass to account for the regular satellites of Jupiter and 
Saturn (only the Saturn results are shown). Furthermore, the ablation and the capture cross-section of meter- 
sized objects can lead to ice/rock fractionation and account for the composition of lapetus, as well as those of 
Titan, Ganymede and Callisto (top right panel). This is the only explanation currently available for lapetus' 
icy composition, and may also indicate that regular satellites overall may be depleted in rock with respect 
to solar composition mixtures, which might explain why their densities are lower than that of similar size 
Kuiper belt objects. 

We stress that while full ice/rock separation is unlikely, only partial separation is needed to account for 
lapetus' icy composition. Such a population can result from the fragmentation of a partially differentiated, 
^^Al-heated first generation of 10 km planetesimals (Mosqueira, Estrada & Charnoz, 2009). Indeed, it is 
likely that most of the mass in the first generation of planetesimals resided in objects 10—100 km in the first 
10^ — 10^ years, so that these objects may have incorporated significant amounts of ^^Al. If so, ~ 10 — 100 
km planetesimals can differentiate (Prialnik & Podolak, 1995) (Phoebe itself may be at least partially dif- 
ferentiated, see Johnson et al. (2009)). A study by Merk & Prialnik (2006) that includes radioactive heating, 
accretion, transformation of amorphous to crystalline ice and melting of water ice in the formation of trans- 
neptunian objects (which take longer to grow than planetesimals forming closer in) with sizes ^ 10 km finds 
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Fig. 1 Left upper panel: the solid-to-gas ratio (dotted line) and surface density of solids (solid line) as a function of radial 
distance from Saturn. Right upper panel: the ice to rock ratio as a function of radial distance. The empty dots correspond to 
Titan and lapetus. The ice mass ratio of solar nebula planetesimal fragments crossing the disk is taken to be 0.3. Left lower 
panel: the mass delivered within a given radius from Saturn. The lower right panel: the planetesimal size that gets either captured 
or ablated as function of radial distance. 
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that the occurrence of liquid water may be common*. The disruption of partially differentiated planetesimals 
would then lead to a population of icy/rocky fragments available to ablate through the extended Kronian gas 
disk. 



2 Alien moons: the capture of Triton and the irregular satellites 

The origin and the subsequent evolution of the irregular satellites have profound implications for our under- 
standing of the satellite systems of the giant planets. In this section we discuss the origins of Triton and of 
the irregular satellites. A separate chapter in this volume (Schubert et al., 2009) extends the discussion to the 
secular evolution of irregular satellites, with a focus on Phoebe. 



2.1 The irregular satellites 

Since the time of their discovery, irregular satellites have been known to belong to a separate population of 
objects. While the regular satellites orbit their parent planets on inner, equatorial and almost circular orbits, 
the orbits of the irregular satellites are about an order of magnitude larger, and are characterised by high ec- 
centricities and inclinations. Moreover, a significant fraction of irregular satelUtes move on retrograde orbits 
(i.e., with inclinations in the range [90° — 180°]). These dynamical features are incompatible with formation 
in prograde circumplanetary disks, which strongly suggests that irregular satellites are captured objects that 
originated elsewhere in the Solar System. 

Depending on the formation regions of the parent bodies, irregular satelUtes would be characterised by dif- 
ferent compositions. Therefore, once captured, irregular satellites would introduce exogenous materials into 
their host systems. Moreover, both theoretical studies and observational results suggest that mass transfer 
between regular and irregular satellites takes place, introducing contaminants on the surfaces of the iimer 
moons. The identification of the irregular satellite capture mechanism (or mechanisms), and of the formation 
regions of the parent bodies, are thus key to our understanding of such contamination processes. 

2.1.1 Capture mechanisms and formation regions 

To permanently capture the parent bodies of the irregular satellites within the Hill radii of the giant planets, a 
fraction of their kinetic energy must be removed during passage. Historically, three dissipation mechanisms 
have been proposed to explain the capture of the irregular satellites: gas-drag (Pollack, Bums & Tauber, 
1979), Pull-Down (Heppenheimer & Porco, 1977), and collisional capture (Colombo & Franklin, 1971). 
The first two scenarios are based on the presence of nebular gas and thus place the capture of the irregular 
satellites at the time of the formation of the giant planets. In the Pollack, Bums & Tauber (1979) scenario gas 
is needed to slow the bodies crossing the Hill's spheres of the giant planets through friction. In the Pull-Down 
scenario by Heppenheimer & Porco (1977) gas is needed to rapidly increase the mass of the giant planets 
and expand their Hill's spheres, thus trapping the irregular satellites in their enlarged gravitational fields. 

* But note that the presence of liquid water is not required for partial differentiation to take place. For instance, enhanced 
conductivity by water vapor flow through porous media may discourage melting, but would still result in a layered internal state 
as the core would become water-depleted (Prialnik et al., 2008). 
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Heppenheimer & Porco (1977) suggest that the gas-drag and Pull-Down mechanisms can be combined into 
a single scenario, where the interplay of the different processes can enhance the capture efficiency. We refer 
the reader to recent reviews (see e.g. Jewitt & Haghighipour (2007) and references therein) for more detailed 
discussions of the viability of gas-based scenarios; nevertheless, we point out that, while such scenarios can 
be appHed to gas-rich Jupiter and Saturn, their applicability to Uranus and Neptune is uncertain. The colli- 
sional capture scenario by Colombo & Franklin (1971) has been somewhat neglected, yet it poses the least 
amount of constraints on the time of capture, has in principle the same capture efficiency when it comes to 
the gaseous and ice/rock giant planets, and may provide an explanation for collisional families among the 
irregular satellites of Jupiter. The Pull-Down scenario by Heppenheimer & Porco (1977) explains the exis- 
tence of these families by allowing for the post-capture disruption of more massive parent bodies through 
collisional events as postulated by Colombo & Franklin (1971). The gas-drag model of Pollack, Burns & 
Tauber (1979) accounts for the existence of collisional families through the break-up of captured bodies due 
to friction. However, owing to the size-dependence of gas drag, one might expect smaller bodies on inner 
orbits and larger ones on outer orbits, which is not observed in the present satellite systems. 
More recently, a new class of dynamical models has been proposed to explain the capture of the irregular 
satellites. These new scenarios fit within the so-called Nice Model, wherein it is argued that the structure 
of the outer Solar System is the result of a phase of chaotic rearrangement of the orbits of the giant planets 
moderated by the presence of a residual disk of planetesimals (for details we refer the reader to the original 
series of articles by Gomes et al. (2005); Morbidelli et al. (2005); Tsiganis et al. (2005)). Since the chaotic 
rearrangement of the orbits of the giant planets would destroy any pre-existing system of irregular satellites 
(see Nesvorny, Vokrouhlicky & Morbidelli (2007), and references within), gas-based models would not be 
viable, with Jupiter as the only possible exception due to its Umited participation in the orbital rearrange- 
ment. Proposed mechanisms to resupply the populations of irregular satellites rely on exchange reactions 
(i.e., the disruption of a binary system of planetesimals while crossing the Hill's sphere of a giant planet, see 
Agnor & Hamilton (2006); VokrouhUcky, Nesvorny & Levison (2008)), and three-body effects due to the 
mutual gravitational perturbations of the giant planets during close encounters in a Nice Model-like scenario 
(Nesvorny, Vokrouhlicky & Morbidelli , 2007). In addition, collisional capture has been studied in more 
detail for Saturn's system in light of the observational data suppUed by Cassini-Huygens mission (Turrini, 
Marzari & Tosi, 2009). 

A detailed discussion of capture mechanisms is beyond the scope of this chapter. In the framework of the 
processes shaping the moons of the giant planets, the most important influence of the irregular satelUtes is 
that of contaminating the surfaces of the inner satellites, introducing exogenous elements with potentially 
different compositional features. The nature of the contaminants delivered depends on the composition of 
the irregular satelUtes, which is Unked to the formation regions of the parent bodies. The underlying as- 
sumption of early studies is that the source of the parent bodies of the irregular satellites is located near the 
orbital regions of the host planets. This would imply that one should expect a compositional gradient with 
heUocentric distance between the families of irregular satelUtes of the different planets, and that the overall 
composition of the irregular satellites should be compatible with that of the regular ones. To the contrary, 
the investigation of the collisional capture scenario performed by Turrini, Marzari & Tosi (2009) shows that 
the possible orbits of the putative parent bodies of the irregular satellites span over a wide range of semima- 
jor axes and eccentricity values, covering the whole outer Solar System. This means that irregular satellites 
orbiting a given planet could have formed in different regions of the Solar System, and thus have significant 
differences in their compositions. The same conclusion also holds true in the Nice Model, due to the orbital 
rearrangement suffered by the early Solar System. 
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2.1.2 Implications of the capture origins 

While our data on the irregular satellites of Uranus and Neptune are still limited, both dynamical and ob- 
servational evidence suggests that collisions between members of the populations of irregular satellites of 
Jupiter and Saturn played an important role in the past. The orbital structure of the Jovian irregular satellites 
bears a clear imprint of the existence of coUisional families (Nesvorny et al., 2003; Nesvorny, Beauge & 
Dones, 2004). The orbital structure of the Kronian irregular satelUtes is more difficult to interpret, due to 
Jupiter's strong gravitational perturbations, yet there are clues of the possible existence of collisional families 
(Turrini, Marzari & Beust, 2008). Moreover, there are strong dynamical indications that Phoebe underwent 
a period of intense collisional activity in the past with now-extinct irregular satellites located near its orbital 
region (Turrini, Marzari & Beust, 2008), and such a hypothesis is supported by the observations of the in- 
struments on-board Cassini of Phoebe's deeply cratered surface (Porco et al., 2005). The colhsional activity 
of the irregular satellites would produce fragments, which could be ejected on orbits that intersect those of 
the regular satellites. Depending on their size, dust grains would be expelled into heliocentric orbits due to 
radiation pressure, or migrate inward due to Poynting-Robertson drag. Dust migrating inward would inter- 
sect the orbits of the regular satelUtes and, depending on their sweeping efficiency, reach different depths in 
the satellite systems. Not all dust grains are collected with the same efficiency: bigger grains migrate more 
slowly, enhancing the probability for the inner satellites to capture them. Also, dust particles in retrograde 
orbits would experience more close encounters with the satellites, and are more likely to be collected than 
grains in prograde orbits. Since the parent bodies of the irregular satellites originated in very different en- 
vironments, this dust production/transfer leads to the contamination of the surfaces of the icy satellites with 
exogenous material characterised by a different chemical composition than the local material, and possibly 
coming from different sources. The origin of the dark material coating the leading hemisphere of lapetus is 
probably the best example of this process. Still, whether the cause of the leading-trailing side dichotomy on 
lapetus is endogenous, as opposed to exogenous, continues to be debated. 



2.2 The case of Triton 

Another case in point is Triton, where the transfer process acts in the opposite way, i.e., the regular satellites 
are the source of contamination on Triton. Like other irregular satellites, Triton's retrograde orbit points to 
a capture origin (Goldreich et al., 1989; Agnor & Hamilton, 2006). Yet, Triton's semimajor axis is over an 
order of magnitude smaller than that of the other irregular satellites, and Triton's size makes it one of the 
most massive moons in the Solar System. The small orbital distance of the satellite from Neptune suggests 
that its close-in, circular orbit is the result of tidal dissipation (Goldreich et al., 1989); it has been pointed 
out, however, that Triton's high inclination {i f« 157°) implies a time-scale of the order of the age of the 
Solar System for tidal dissipation to shrink the orbit of the satellite from probable post-capture values of 
the semimajor axis to its present location (Cuk & Gladman, 2005). Triton's gravitational perturbations from 
its initial post-capture orbit would have caused pre-existing regular satellites to collide with each other and 
shatter. Taking the mass of the pre-existing satellites to be comparable to that of Triton, which these authors 
justify by analogy with the Uranian satellite system, the resulting debris disk could have shrunk Triton's orbit 
through scatterings and colhsional debris-drag on a time-scale of about 10^ years, possibly allowing for the 
survival of the outermost regular and irregular satellites. Such a hypothesis might be in agreement with the 
results obtained by Schenk & Zahnle (2007) in reprocessing the Voyager's data on Triton, suggesting that 
the asynmietric distribution of craters on the surface of the satelUte could be best explained by assuming 
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a planetocentric source for the impactors, i.e., that the craters on Triton's leading hemisphere are due to 
impacts with a population of prograde small satellites or collisional shards. The paucity of craters on Triton 
could be an indication that the surface of the satelhte is extremely young, with its age varying in the range 
10^ — 10^ years, diagnostic of an active geology possibly Unked to the tidal evolution of the satellite. This 
young surface age could indicate that the planetocentric impactors responsible for the craters could be the 
surviving members of the original debris disk. While the possibihty of hnking Triton's origin and evolution 
to its present state is intriguing, it should be noted that at the moment no quantitative study of the collisional 
lifetime of such a debris disk exists. Moreover, the data available on Triton and on Neptune's system in 
general are quite limited, making it difficult to derive better observational constrains for the theoretical 
models. Nevertheless, it is clear that one needs to take into account Triton's history and its exogenous origin 
to properly interpret the data currently available and those to be supplied by future space missions. While 
the Voyager's data indicate that Triton's surface is quite young and possibly uncontaminated by external 
material, the craters on the leading hemisphere could contain residual material from the impactors, perhaps 
more similar in composition to the present regular satellites of Neptune than to the uncontaminated surface 
of Triton. At any rate, Triton's capture disrupted Neptune's pre-existing satellite system, implying that care 
needs to be exercised when attempting to extract information about proto-Neptune's subnebula from the 
structure of the present satellite system of the giant planet. 



3 Evolution: the path to present 

As planetary formation sets the boundary conditions for satellite formation, so satellite formation sets the 
initial conditions for satellite evolution. Dynamical and geochemical evolution, in fact, secularly reshapes the 
satellites, erasing or modifying their primordial features. Therefore, in order to enhance our understanding 
of the present state of the moons of the outer Solar System, we need to be able to discriminate the effects of 
their secular evolution from those properties inherited from their time of formation. Nevertheless, thanks to 
the recent Cassini mission to Saturn and the efforts of those who made it possible, we currently have in hand 
a sufficient number of reliable observational constraints to give us justifiable hope that it is indeed possible 
to make significant progress in our understanding of the origin and evolution of the satellite systems of the 
giant planets. 
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